Self-Interacting Dark Matter with Flavor Mixing 
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The crisis of both cold and collisional dark matter (DM) models is resolved by postulating flavor 
mixing of DM particles. Flavor-mixed particles segregate in the gravitational field to form dark halos 
composed of heavy mass eigenstates. Since these particles are mixed in the interaction basis, elastic 
collisions convert some of heavy eigenstates into light ones which leave the halo. This annihilation- 
like process will soften dense central cusps of halos. The proposed model accumulates most of the 
attractive features of self-interacting and annihilating DM models, but does not suffer from their 
severe drawbacks. This model is natural does not require fine tuning. 
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Dark matter constitutes most of the mass in the Uni- 
verse, but its nature and properties remain largely un- 
known, except, perhaps, it is cold (at the time of decou- 
pling, i.e., massive) and dark (weakly interacting). Re- 
cent observations indicate a failure of the conventional 
cold dark matter (CDM) paradigm on small (galactic) 
scales. The resolution of the discrepancy between the- 
ory and observations calls for new ideas and should allow 
to further constrain properties and nature of otherwise 
undetectable dark matter. 

A model of structure formation in a universe with 
CDM is in excellent agreement with observations on large 
scales (3> Mpc) 0] which, thus, supports the hypothesis 
that dark matter particles are heavy (> GeV) and weakly 
interacting with baryonic matter and photons. Recent 
observations on galactic and sub-galactic scales appear 
to be in conflict with the CDM model, which predicts 
singular (p ~ 1/r or even steeper) central density cusps 
[| H and a large number of low- mass halos || |). 
In contrast, observations indicate low-density cores with 
shallow profiles || 0, || and a smaller number of satel- 
lites H- Also, some observations seem to indicate the 
invariance of the core density with respect to the mass of 
the halo |fto). 

A number of models have recently been proposed in 
order to explain the discrepancy. Some of them, such as 
the models of self-interacting, warm, annihilating, and 
decaying dark matter 12, [l4| are just (slight) 



modifications of the standard CDM paradigm, while oth- 
ers, e.g., the fluid, repulsive, fuzzy, bosonic, mirror, etc. 
dark matter models Jl|, [l§ [H], [l§ fl9j are rather exotic. 
Nearly all theories may be tuned to solve the problem, 
they may be motivated from particle physics, and, finally, 
they cannot be ruled out by present observational data. 

The fact that many predictions of the CDM cosmolog- 
ical models are in excellent agreement with a wide va- 
riety of observational results, from the microwave back- 
ground fluctuations to the large-scale structure, suggests 
that the CDM model is a good first approximation, but 
it has to be corrected on a small (galactic) scale. Among 
many possibilities, the simplest one is to assume that 
weakly interacting CDM particles are, in fact, strongly 



self- interacting O, |TJ . Evolution of the structure in 
the presence of self-interacting dark matter (SIDM) in- 
deed results in a smaller number of satellites and forma- 
tion of constant density cores within a dynamical time- 
scale modified by collisions |23, 



26 . On a 



24 



longer time-scale, however, a SIDM halo is unstable to 
core collapse, which leads to the formation of a steep 
~ 1/r 2 density cusp, making the problem even worse p5[ ] 
(but see Ref. ^6|). The detection of elongated cores in 
clusters via lensing observations p7| is inconsistent with 
the simple SIDM model. Thus, it seems that such an 
one-component model is an over-simplification, because, 
in general, dark matter particles may be flavor mixed, 
that is their interaction (or flavor) eigenstates are super- 
positions of mass eigenstates. 

In this paper we study the properties of dark matter 
halos in the model of the SIDM with flavor mixing. We 
demonstrate that elastic scattering of such dark matter 
particles results in a conversion of mass eigenstates and 
calculate the cross-section of this process. The conver- 
sion process will lower the density in the core region and 
prevent the cusp formation. This is somewhat similar to 
the annihilating dark matter model [ jl3| , but our model 
does not suffer from the "annihilation catastrophe" in 
the early Universe because of the dynamical equilibrium 
of the opposite conversion processes in a homogeneous 
medium. 

For the sake of simplicity, let us assume that there ex- 
ist two flavors, the strongly self- interacting, |2), and the 
non-self-interacting or "sterile", |<S), ones. These inter- 
action eigenstates are a superposition of two mass eigen- 
states, the "heavy", \rrih), and the "light", \mi), ones. 
We write 
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where i? is the mixing angle. Throughout the paper, we 
assume for simplicity that m/j 3> mj , so that heavy states 
are non-relativistic and light states are relativistic. Gen- 
eralizations of the above simplifications are straightfor- 
ward. We should note that, in general, the treatment 
of non-relativistic species requires the wave packet ap- 
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proach [£8j, an d not the plane wave one used above. In 
some instances, this may result in measurable corrections 
p9| , p0[ . Here we neglect by these complications and use 
the standard plane wave approach, keeping in mind that 
a particle is a localized in space wave packet with the 
width Sx determined by the production process. 

The concept of flavor eigenstates arises when interac- 
tions of particles are considered. In the field-free theory 
the particle fields in the mass basis have a physical mean- 
ing instead. In general, such mass eigenstates have differ- 



ent momenta and energies, Ef il = \p h 
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tween the interactions, the mass eigenstates propagate in- 
dependently, with different velocities \h,i = Ph,ic 2 / En,i- 
Thus at times 
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these states are separated from each other and their wave 
functions no longer overlap. The separation time above 
is negligibly small compared to a galactic dynamical time 
scale; thus, this process is essentially instantaneous. In a 
gravitational field different eigenstates segregate by mass. 
Non-relativistic |m^,}-states form halos and the large- 
scale structure, they represent the conventional cold dark 
matter. Relativistic |m;)-states correspond to the "hot" 
dark matter component which cannot form structures. 
This component has a much smaller total mass, however, 
(the number density of both mass states are equal by 
symmetry) and, hence, is dynamically unimportant. 

As dark matter halos form, the density in the central 
parts increases and, at some point, self-interactions of the 
dark matter particles become important. Let us consider 
the elementary act of elastic scattering of two particles, 
i.e., \rrih) eigenstates. The initial wave function of two in- 
teracting particles in the center of mass frame, according 
to equation ([!]), is 
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where "+" corresponds to ^ symmetric to interchange of 
particles (integer total spin) and "— " - to an antisymmet- 
ric (half- integer total spin), the exponents represent 
two waves, propagating to the right, and to the left, and 
e ±ikz _ ^ e ikz _|_ e -jfcz^ j s short-hand notation. For 

scattering, the interaction basis is appropriate, rather 
than the mass basis, hence the expansion above. Dur- 
ing the scattering event, only (X)-component is chang- 
ing, since |iS)-component does not interact. The wave 
function at large distances after scattering thus becomes 
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where the combination f± (0) — f(6)±f{n — 0) arises be- 
cause particles are indistinguishable and f(8) is the am- 
plitude of scattering of flavor states. The radial part of 
the wave function represents a diverging scattered wave. 
One can clearly see that an initial heavy eigenstate ac- 
quires, upon scattering, a light eigenstate admixture. In 
other words, a heavy eigenstate may be converted into 
a light eigenstate. The differential cross-section of the 
conversion process is the scattering amplitude squared, 
that is 



dcr c ,± = Icosfisiwd f±(9)\ dQ, 



(5a) 



where dfl = 2irsm9dd is the solid angle. The cross- 
section of pure scattering is 



da 



s,± 



|cos 2 tf/±((9)| dfl. 
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Of course, energy and momentum must be conserved in 
these processes. Introducing the integral cross-section as 
a± — J \ f±(9)\ 2 dtl, we write 



er c .± = sin -d cos -d a±, 



<j s ,± = cos "d a±. 



(G) 



It is assumed here that the interacting particles of spin s 
have a particular value of their total spin. If it is not so, 
one has to average over all spin states, which yields |3l[] 
for the conversion cross-section: 



f 
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and similarly for the scattering cross-section cr s . The 
cross-sections for the light mass eigenstate interactions 
may be obtained similarly. 

Let us now consider a dark matter halo. When dark 
matter particles (heavy mass eigenstates) interact with 
each other, they can either scatter with some cross- 
section <j s or be converted into light mass eigenstates 
with the cross-section 



a c = tan ■& <r s 



(8) 



Since these light eigenstates are relativistic, the conver- 
sion process is analogous to the annihilation of dark mat- 
ter particles. Since the emergent light eigenstates are 
relativistic, they leave the dark halo, thus, decreasing its 
density and mass. If the mixing angle is not too small, 
■d ~ I, collisional heat transport, which flattens the core 
region by making it nearly isothermal, and dark matter 
"annihilation" (conversion) become important at roughly 
the same density of particles and, hence, operate simul- 
taneously. Hence, on the onset of the core collapse, when 
the infall velocity is small, there is enough time to an- 
nihilate the density enhancement and prevent formation 
of a central cusp. Only direct particle simulations are 
capable of determining what the spatial structure of the 
dark matter halo will be in this regime. 
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As one can see, most of the assumptions adopted in the 
beginning are not crucial for our model. The model re- 
quires only that the interaction strengths of flavor eigen- 
states and the masses of mass eigenstates be (signifi- 
cantly) different. No fine tuning is necessary. 

A remarkable difference of the eigenstate conversion 
process from pure annihilation is that particles are not 
permanently destroyed via conversion. This is a re- 
versible process: \rrih} eigenstates are "regenerated" in 
interactions of \mi) eigenstates. By symmetry, the con- 
version cross-sections of regeneration and annihilation 
are identical. Now consider the early Universe. The 
density fluctuations in it are negligible. Therefore, both 
processes occur at the same rate. They balance each 
other exactly, no matter how large the particle density 
is. Thus, no excess of particles of one type over particles 
of another type is produced. In contrast, conventional 
annihilation is extremely rapid at high densities. The an- 
nihilating dark matter model [ |l3| , thus, suffers from the 
"annihilation catastrophe" . Namely, dark matter must 
be destroyed in the early Universe epoch, long before the 
large scale structure began to form. 

We also should note that the conversion and pure an- 
nihilation cross-sections have, in general, different depen- 
dence on the relative velocity, v, of interacting parti- 
cles. In the absence of long-range interactions and in 
the non-relativistic limit, the annihilation cross-section 
scales usually as a a oc v~ , while cr c , being proportional 
to the cross-section of elastic scattering, does not depend 
on v, i.e., a c oc v°. These scalings are typical but not 
universal, however, and depend on the interaction model 
at hand @ || ||. 

The dynamics of SIDM model with flavor mixing is 
quite rich. In the limit of no mixing, $ = 0, the model 
reduces to the standard SIDM model. If the mixing angle 
is rather small, elastic collisions rapidly establish isother- 
mal equilibrium throughout a dark matter halo. The halo 
becomes relatively spherical in this case. Moreover, its 
core will collapse. Particle conversions will be important 
deep inside the collapsed core where the density is high 
enough. Finally, if 1} > 1, conversions of dark matter par- 
ticle rapidly lower the halo density. An almost spherical 
isothermal core never forms. The detection of elongated 
shapes from lensing observations of clusters |^7j seems to 
favor the last scenario. 

To summarize, we presented a model of self-interacting 
dark matter with flavor mixing. This model is well moti- 
vated from particle physics (e.g., mixing of quarks, neu- 
trinos, etc.) within the standard model, it is robust and 
does not require fine tuning of any kind. We have demon- 
strated that dark matter halos are composed of particle 
mass eigenstates. Being mixed in the interaction basis, 
some of heavy eigenstates are converted in elastic colli- 
sions into light eigenstates and leave the halo, in a strik- 
ing analogy with halo annihilation. Accumulating many 
attractive features of SIDM and annihilating dark matter 



models (e.g., formation of soft cores and a fewer number 
of sub-halos, while preserving CDM properties on a large 
scale), our model does not suffer from their caveats, such 
as the core collapse induced r _2 -cusps and too spherical 
cores in the SIDM model, and the catastrophically strong 
annihilation of dark matter in the dense early Universe. 

The author is grateful to Jim Peebles for interesting 
and insightful discussions and his interest in this work 
and to CITA for the financial support. 
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